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Abstract Photometric analysis of BVRc light curves 
of newly discovered eclipsing binary GSC 0008-00901 is 
presented. The orbital period is improved to 0.28948(11) 
days. Photometric parameters are determined, as well. 
The analysis yielded to conclusion that system is an 
over-contact binary of W UMa type with components 
not in thermal contact. The light curves from 2005 show 
the presence of a spot on the surface of one of the com- 
ponents, while light curves from 2006 are not affected 
by maculation. 
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1 Introduction 

The variability of CCS 00008-00901 (a 20 oo= 00 h 13 m 
22.689 s , <$ 2000 =+05° 40' 09.25"), located in the field 
of eclipsing binary DV Psc, was recently reported by 
Parimucha (2007). He analyzed V and Rc light curves 
and showed that the system is most probably contact 
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binary with spotted regions on one of the components. 
Cutri et al. (2003) reported this star as suspected for 
variability during their 2MASS project. 



2 Observations and data reduction 

Our photometric observations of the DV Psc field has 
started on September 4, 2005 and ran to October 26, 
2006. The system was detected on CCD images from 5 
nights in 2005. Observation field in 2006 was selected 
in the manner to obtain system during all observations 
of DV Psc. All these observations were obtained at the 
Stara Lesna (SL) Observatory of Astronomical Institute 
of the Slovak Academy of Sciences and at the Kolonica 
Observatory (KO). The journal of CCD observations is 
given in Tab. 1. 

For observations from SL observatory the 50cm tele- 
scope equipped with the SBIG ST-10 MXE CCD cam- 
era with standard BVRc filters was used. Observations 
from KO were performed by 265/1360 mm Newton type 
telescope with Meade DSI Pro CCD camera with Sony's 
ExView HAD Monochrome CCD Image Sensor. No fil- 
ter was used for these observations. 

The CCD images were reduced in the usual way 
(bias and dark subtraction, flat-field correction) in soft- 
ware package SPHOTOM developed by the first author. 
This package was also used for aperture photometry 
of all detected stars on the CCD frames. This pro- 
cedure is based on the latest SExtrator code (Bcrtin 
and Arnouts, 1996). The brightness of the star was 
determined with respect to BD +04°19 (GSC 00008- 
00949) with B = 12.04 mag, V = 11.47 mag, R = 11.22 
mag, (Monet et al., 2003). The stability of comparison 
star was tested with respect to GSC 00008-00743 and 
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Table 1 The journal of CCD observations of CSC 0008-00901 
obtained at Stara Lesna (SL) and Kolonica (KO) observatories. 
The phase interval is determined according to ephemeris (2). Ob- 
servations denoted by * were used only for minima times deter- 
mination. 



Date 


Time (UT) 


Phase 


Filter 


Obs. 


Sep 04, 05 


23:12 


- 03:04 


0.221 


- 0.598 


BVR 


SL 


Sep 23, 05 


19:45 


- 22:03 


0.200 


- 0.527 


BVR 


SL 


Sep 26, 05 


20:39 


- 01:01 


0.696 


- 0.330 


BVR 


SL 


Oct 04, 05 


18:36 


- 01:22 


0.097 


- 0.954 


BVR 


SL 


Oct 27, 05 


19:21 


- 23:26 


0.615 


- 0.197 


BVR 


SL 


Aug 15, 06 


22:35 


- 02:25 


0.769 


- 0.320 




KO* 


Aug 17, 06 


22:28 


- 02:26 


0.661 


- 0.234 




KO* 


Aug 26, 06 


21:28 


- 02:43 


0.608 


- 0.363 




KO* 


Sep 16, 06 


19:02 


- 00:45 


0.290 


- 0.856 




KO* 


Oct 17, 06 


18:04 


- 00:54 


0.751 


- 0.713 


BVR 


SL 


Oct 18, 06 


19:27 


- 00:57 


0.405 


- 0.192 


BVR 


SL 


Oct 26, 06 


18:35 


- 00:12 


0.914 


- 0.714 


VR 


SL* 



e 14. 



R 13.0 




Phase 



Table 2 Weighted average BVR times of the primary (I) and 
secondary (II) minima of CSC 00008-00901. The standard errors 
of the minima are given in parentheses. 



HJD 


Type 


Filter 


2400000+ 






53618.5935(2) 


II 


BVR 


53640.4491(1) 


I 


BVR 


53648.4064(2) 


II 


BVR 


53671.4278(2) 


I 


BVR 


53963.5096(1) 


I 


none 


53965.5402(1) 


I 


none 


53974.5140(2) 


I 


none 


53995.3626(5) 


I 


none 


54026.3318(2) 


I 


BVR 


54026.4709(2) 


II 


BVR 


54027.3427(2) 


II 


BVR 


54027.4921(1) 


I 


BVR 


54035.3076(2) 


I 


VR 


54035.4492(1) 


II 


VR 



BD +04°18 and was found to be stable within -0.003 
mag during all observations. 

Instrumental magnitudes b, v and r (for SL observa- 
tions) and Am (for KO observations) were determined. 
Nightly atmospheric extinction coefficients were com- 
puted from the comparison stars measurements and the 
correction for an atmospheric extinction was applied 
to these observations. Above-atmosphere instrumental 
magnitudes b ,v and r were determined and consec- 
utively were corrected to the standard BVRc interna- 
tional system using following equations : 



V = v + k v (B-V) 
B - V = k bv (ba - vq) 
V - R c = k vr (v - r ), 



(1) 



where k v —0, fc bt) =0.834 and fc vr =1.467. Average trans- 
formation coefficients are routinely determined one or 
twice per month from measurements of the Landolt 
standarts (Landolt, 1992). KO observations were not 
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Fig. 1 The BVR C light curves of GSC 00008-00901 obtained in 
2005 (top) and 2006 (bottom). 



reduced to international system and were used for min- 
ima times determination only (see Section 3). 

GSC 00008-00901 is relatively faint (with V max ~ 
13.9 mag) star and exposure times of our observations 
in different passbands were selected in focus to DV PSc. 
It strongly affected a quality of the obtained light curves 
of this system with our instruments. Because of poor 
quality of photometric data, we had to remove from 
light curves analysis SL observations from October 26, 
2006. This observation was suitable only for minima 
times determination. 

The average accuracy of the observations in B and 
V passband is —0.005 mag and in Rc passband —0.008 
mag. The accuracy of the KO observations is about 0.02 
mag. 
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3 Minima times and period determination 

Our observations enabled us to determine totally 14 
minima times of GSC 00008-00901. They are listed in 
Tab. 2. The times of minima were determined sepa- 
rately for all filters (SL observations) using Kwee and 
van Woerden (1956) method and then weighted aver- 
age minima from all filters were determined. For com- 
putation of minima times we have used only data in 
a phase interval ±0.05 around the minimum. This ap- 
proach minimize the influence of the minima asymme- 
tries (van't Veer, 1973). 

Using period analysis of the light curves during dis- 
covery process, the period P — 0.289437(29) days was 
found. A brief inspection of the light curve obtained on 
October 4, 2005 showed that observed minimum is sec- 
ondary. It enables us to estimate minima types as listed 
in Tab. 2. 

Listed minima times were used to improve period 
and provide linear ephemcris for primary minimum: 

Min I = 2453648.26421 + 0.28948 xE, 

±23 ±11 ^ ' 

which was used to phase all our data. 

The BVR C light curves of GSC 00008-00901 ob- 
tained in 2005 and 2006 seasons, phased with ephemeris 
([2]) are depicted in Fig. 1. 

4 Light curves analysis 

A more detailed look at light curves of GSC 00008- 
00901 obtained in 2005 and 2006 seasons shows that 
they are quite different. Observations from 2006 show 
symmetric light curves in all passbands with brightness 
in the both maxima at the same level. On the other 
hand, data from 2005 indicate that the brightness of the 
system in maximum in phase 0.75 (MaxII) is brighter 
than maximum in phase 0.25 (MaxI).To quantify this 
effect, we performed quadratic polynomial fit of the ob- 
served data around the MaxI and MaxII and we found 
that MaxI-MaxII is 0.033 mag for B filter, 0.028 mag 
for V and 0.02 mag for Rq filter, respectively. This indi- 
cates a presence of a spot(s) on the surface of the stars. 
So, we decided to solve light curves from this seasons 
separately. 

The analysis of the light curves was performed by 
software package PHOEBE (Prsa and Zwitter, 2005). It 
is a modeling package for eclipsing binary stars, which is 
based on widely used WD program (Wilson and Devin- 
ney, 1971; Wilson, 1979; 1990). 

In the process of light curve analysis of the binary we 
have adopted and fixed following parameters: the mean 
temperature of the primary component, T\ = 5350X, 
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Fig. 2 The variations of the normalized weighted sum of squares 
of residuals w(0 — C) 2 for all light curves from 2006 for dif- 
ferent values of mass ratio q and orbital inclination i. 

based on the B — V=0.72±0.09 mag color index, which 
corresponds to G8 spectral type (Cox, 2000), coeffi- 
cients of gravity darkening g\ = gi = 0.32 (Lucy, 1967) 
and the bolometric albedo coefficients of A\ — — 
0.5 (Rucinski, 1969), which are appropriate values for 
the convective envelopes. A Kurucz (1993) model of 
stellar atmospheres was applied to the stars assum- 
ing solar composition. The second order bolometric and 
monochromatic limb darkening coefficients for logarith- 
mic law, were interpolated from van Hamme (1993) ta- 
bles. 

The weights of individual data points were estab- 
lished as 1/er 2 , where a is standard error of point de- 
rived during photometric measurement. In our solu- 
tions, we used all data points, we have not calculated 
normal points as usual. Using this approach we do not 
loose information in our data. 

We fitted all three light curves simultaneously. Based 
on appearance of the light curves, we set PHOEBE so- 
lution to overcontact binary, which is not in thermal 
contact. Initial values of free parameters were deter- 
mined manually to give a visually good fit with our 
data. With these values, a correct solution was quickly 
found. We have run differential correction fit as long as 
output corrections were smaller than the errors of the 
fitted parameters. 

In the first step, we solved light curves from 2006, 
when we assumed that there are no spots on the surface 
of the stars. According to Binncndijk (1957) dividing 
of W UMa systems into A and W types, we concluded 
that this system is a W-type. because of G8 spectral 
type and short orbital period. It means, that primary 
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Table 3 Photometric elements of GSC 00008-00901 and their 
standard errors (in parenthesis) for 2006 light curves: i - orbital 
inclination, q - photometric mass ratio, f! - surface potential, / 
- fill-out factor, ri,r2 -volume mean fractional radii, T\ - fixed 
temperature of the primary component, T2 - temperature of sec- 
ondary component, ^ui(0 — C) 2 - normalized weighted sum of 
squares of residuals for all light curves 



Par? 



Table 4 Temperature of the secondary component T2 and spot 
parameters and their standard errors (in parenthesis) for differ- 
ent configurations: CP - cool spot on primary, CS - cool spot on 
secondary, HP - hot spot on primary, HS - hot spot on secondary, 
Ion - longitude of the spot, lot - latitude of the spot, rad - ra- 
dius of the spot, k - temperature factor (fixed), ^ui(0 — C) 2 
- normalized weighted sum of squares of residuals for all light 
curves 



Parameter 


CP 


CS 


HP 


HS 


T 2 [K] 


5158(8) 


5171(5) 


5162(7) 


5169(9) 


lon[°] 


282(5) 


112(4) 


105(6) 


285(5) 


lat[°] 


71(13) 


73(14) 


70(14) 


71(11) 


rad[°] 


26(2) 


30(3) 


22(3) 


21(2) 


k 


0.95a 


0.95a 


1.05a 


1.05a 


}]w(o-cy 


2.235 


2.214 


2.236 


2.229 



minimum is caused by occultation of the less massive 
and hotter primary star by more massive but cooler 
secondary component. Beacause of W-type definition, 
the mass ratio q > 1. 

To find correct values of mass ratio q and orbital 
inclination i, we computed normalized weighted sum of 
squares of residuals ^ w(0 — C) 2 for fixed values of q in 
interval 1.0-1.8 and i in interval 60-90° and looked for a 
global minimum. A detailed view of the variations of the 
sum of squares around the global minimum is depicted 
at Fig 2. We obtained photometric parameters, which 
are listed in Tab. 3. The best fits with corresponding 
O-C for all passbands are depicted in Fig. 3 (bottom). 

These parameters were used in the next solution of 
2005 light curves. As was mentioned before, the light 
curves from this season are affected by maculation. The 
process of the searching for the spot's parameters is 
complicated by the fact that observed O'Connell ef- 
fect could be caused by different combinations of the 
spot's location (on primary or secondary component) 





72.35(15) 


0.8 - 


q[°] 


1.332(11) 


0.7 - 


Q 


4.194(15) 


0.6 - 


f 


0.114(26) 




n 


0.359(5) 




r 2 


0.417(9) 




Ti[K] 


5350a 




T2IK} 


4995(8) 




L 1 /(L 1 +L 2 ) B 


0.557(30) 




Li/(Li + L 2 ) V 


0.530(28) 




Li/(Li +L 2 )r c 


0.513(25) 




>>(0-C) 2 


0.859 
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Fig. 3 The best fits of BVR C light curves of GSC 00008-00901 
for 2005 data (top) and 2006 data (bottom). B and V data are 
shifted for a clarity. 



and spot's temperature (hot or cool spot). To quan- 
tify this effect, we tried to fit light curves from 2005 
for 4 configurations of the spot location and tempera- 
ture as listed in Tab. 4. In these solutions we have fixed 
temperature factor (ratio between temperature of the 
spot and surrounding atmosphere) of the spot as well 
as parameters from 2006 solution except temperature 
of secondary component. As could be seen from values 
of J2 W (0 ~ C) 2 in Tab. 4, to select correct solution 
is practically impossible. In Fig. 3 (top) we show light 
curve solution in the case when a cool spot is located on 
the secondary component. This solution has a minimal 
weighted sum of squares of residuals for all light curves. 
The corresponding 3D model of the system with a cool 
spot on secondary component is shown in Fig. 4. 
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Fig. 4 3D model of GSC 00008-00901 with a cool spot on the 
secondary component 



5 Discussion and Conclusion 
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Our light curve analysis of BVRc photometric observa- 
tions of recently discovered eclipsing binary GSC 00008- 
00901, which is located in the field of another eclipsing 
binary DV Psc showed, that this system is an overcon- 
tact binary of W UMa type. Considering of its spectral 
type and orbital period, we concluded that the primary 
minimum is occultation (primary component is eclipsed 
by secondary one). So, this system belongs to an W- 
type of W UMa binaries. The temperature difference 
between components is about 350K. 

In spite of problems with quality of the observa- 
tions we showed that light curves from the both seasons 

2005 and 2006 are different. The light curves from 2005 
show different brightness in the both maxima, while 

2006 light curves arc symmetric. Asymmetry of the light 
curves from 2005 could be explained by the spot on 
one of the components, however, some conclusions on 
spot's temperature and/or its location (on primary or 
secondary component) are uncertain. 
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